Colleen B. Bove1,2*, Laura Mudge1, and John F. Bruno1

1 The Department of Biology, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, 27599-3280 USA.

2 The Department of Biology, Boston University, Boston, Massachusetts, 02215 USA.


Abstract:

Abstract text will go here


Citation:

Citation will go here


Caribbean Ecoregion Coral Reefs

Figure 1

## OGR data source with driver: ESRI Shapefile 
## Source: "/Users/colleen/Dropbox/Git/CaribbeanSST/data/EcoRegions", layer: "Caribbean_ecoregions"
## with 17 features
## It has 9 fields

Caribbean coral reef site locations and ecoregion designation. The colour of each reef represents the designated ecoregion (Spalding et al. 2010) and n denotes the number of unique reefs within that ecoregion.


SST Extraction and Analysis

Figure 2

Historic SST trend on Caribbean coral reefs (1871–2020). Long-term SST records (HadISST) on Caribbean coral reefs depicting A) mean monthly SST each year (represented by line colour: blue to red) and B) GAM smoothed annual mean SST time (black line), annual maximum (red line), and annual minimum (blue line) SST. The grey dashed horizontal line denotes the overall mean SST for all sites over the entire period (27.3 °C) and the grey ribbon represents the 95% confidence interval around the true annual SST mean through time.


Figure 3

Historic SST trends on coral reefs within ecoregions (1871–2019). Long-term SST records (HadISST) on Caribbean coral reefs separated by ecoregion depicting GAM smoothed annual mean SST time (black line), annual maximum (red line), and annual minimum (blue line) SST. The grey dashed horizontal line denotes the mean SST over the entire period and the grey ribbon represents the 95% confidence interval around the true annual SST mean for each ecoregion.


Figure S2

Historic SST records (1981-2019; Pathfinder) on Caribbean coral reefs depicting A) mean monthly SST each year (represented by line colour: blue to red) and B) GAM smoothed annual mean SST time (black line), annual maximum (red line), and annual minimum (blue line) SST. The grey dashed horizontal line denotes the overall mean SST for all sites over the entire period (27.23 °C) and the grey ribbon represents the 95% confidence interval around the true annual SST mean through time.


Figure S3

Comparison of HadISST (1871-2020) and Pathfinder (1981-2019) SST recorded on Caribbean coral reef locations. The high-resolution Pathfinder is represented as darker data over the long-term HadISST. Both datasets are represented by GAM smoothed annual mean SST time (solid line), annual maximum (red line), and annual minimum (blue line) SST. The dashed horizontal line denotes the overall mean SST for all sites over the entire period and the grey ribbon represents the 95% confidence interval around the true annual SST mean through time.


Figure S4

Still need to update this figure


HadISST GIF

Long-term SST records (HadISST) on Caribbean coral reefs depicting mean monthly SST each year (represented by line colour: blue to red).


Caribbean Basin Maps and Marine Heatwaves

Figure 4

Warming patterns throughout the Caribbean Sea. Increasing warming events across the Caribbean depicted through A) rate of SST change (°C per decade) from 1981 to 2019 (Pathfinder; mean slope 0.23 ± 0.087 °C per decade) and B) increasing marine heatwave events (slope of counts per year). Grey ocean area was not included in these analyses.



Figure 5

Need to add MHW metrics here



Figure S5

Significance of A) rate of SST change (°C per decade) and B) number of marine heatwave events per year across the Caribbean depicted in Figure 4. Grey ocean area was not included in these analyses.



Figure S6

A) Rate of SST change (°C per decade) over the duration of the HadISST database across the Caribbean from 1871 to 2020 (mean slope 0.04 ± 0.014 °C per decade) and B) significance of rate of SST change. Grey ocean area was not included in these analyses.



Figure S7

A) Rate of SST change (°C per decade) across the Caribbean from 1981 to 2020 (HadISST; mean slope 0.16 ± 0.054 °C per decade) and B) significance of rate of SST change. Grey ocean area was not included in these analyses.



Supplemental Tables

Table S1

Table S1. Marine heatwave properties examined in this study

Metric Units Description
Frequency Number of events Number of discrete mhw events
Total MHW days days The total number of days a location experienced a MHW per year
Duration days The number of days between the start and end date of each distinct MHW event
Peak intensity °C The maximum temperature, above the seasonal varying climatological mean, reached during the MHW event
Onset rate °C/day “Rate of temperature change between onset date and date of peak intensity”
Return time days The number of days elapsed since a previous mhw event in that location

Table S2

Table S2. Calculated warming rates from both HadISST and Pathfinder databases for different temporal ranges.

Temperature parameter HadISST (1871-2020) HadISST (1990-2020) Pathfinder (1981-2019) Pathfinder (1990-2019)
Caribbean Basin (°C per decade) 0.04 0.16 0.23 NA
Caribbean Basin (total °C for period) 0.60 0.50 0.90 NA
Caribbean Reefs (°C per decade) 0.04 0.17 0.19 0.2
Caribbean Reefs (total °C for period) 0.60 0.53 0.72 0.62

Table S3

Table S3. Mean ocean warming rate (°C per decade; with 95% confidence interval) and total increase in temperature (°C) on coral reefs within each Caribbean ecoregion since the noted year of inflection point. The inflection point for each ecoregion was identified as the year in which annual warming rates noticeably increased based on the GAM smoothed annual means (see Figures 3 and S4). Years indicated with an asterisk (*) represent inflection points identified for ecoregions that occurred before the Pathfinder dataset record, therefore those rates and total warming values were calculated from the beginning of the record in 1981.

Ecoregion Inflection point Mean rate (°C per decade) Lower confidence interval Upper confidence interval Warming extent (°C)
Bahamian 1986 0.16 0.151 0.170 0.53
Eastern Caribbean 1985 0.26 0.255 0.262 0.88
Floridian 1981 0.20 0.190 0.204 0.76
Greater Antilles 1986 0.16 0.160 0.168 0.53
Gulf of Mexico 1981 0.21 0.186 0.241 0.80
Southern Caribbean 1981 0.26 0.258 0.268 0.99
Southwestern Caribbean 1981 0.19 0.182 0.192 0.72
Western Caribbean 1993 0.13 0.126 0.138 0.34

Session information

All code was written by Colleen B. Bove, feel free to contact with questions.

Session information from the last run date on 2021-03-16:

## R version 3.6.3 (2020-02-29)
## Platform: x86_64-apple-darwin15.6.0 (64-bit)
## Running under: macOS Catalina 10.15.7
## 
## Matrix products: default
## BLAS:   /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRlapack.dylib
## 
## locale:
## [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
## 
## attached base packages:
## [1] grid      stats     graphics  grDevices utils     datasets  methods  
## [8] base     
## 
## other attached packages:
##  [1] repmis_0.5          kableExtra_1.3.1    plotly_4.9.3       
##  [4] foreach_1.5.1       mgcv_1.8-33         nlme_3.1-151       
##  [7] gganimate_1.0.7     fields_11.6         spam_2.6-0         
## [10] dotCall64_1.0-0     binr_1.1            RColorBrewer_1.1-2 
## [13] xts_0.12.1          zoo_1.8-8           rnaturalearth_0.1.0
## [16] forcats_0.5.0       stringr_1.4.0       dplyr_1.0.2        
## [19] purrr_0.3.4         readr_1.4.0         tidyr_1.1.2        
## [22] tibble_3.0.4        tidyverse_1.3.0     cowplot_1.1.1      
## [25] rgdal_1.5-19        maptools_1.0-2      viridis_0.5.1      
## [28] viridisLite_0.3.0   raster_3.4-5        sp_1.4-5           
## [31] ggrepel_0.9.0       ggplot2_3.3.3       ncdf4_1.17         
## [34] sf_0.9-7           
## 
## loaded via a namespace (and not attached):
##  [1] colorspace_2.0-0        ellipsis_0.3.1          class_7.3-17           
##  [4] fs_1.5.0                rstudioapi_0.13         farver_2.0.3           
##  [7] fansi_0.4.1             lubridate_1.7.9.2       xml2_1.3.2             
## [10] codetools_0.2-18        splines_3.6.3           R.methodsS3_1.8.1      
## [13] knitr_1.30              jsonlite_1.7.2          broom_0.7.3            
## [16] dbplyr_2.0.0            R.oo_1.24.0             rgeos_0.5-5            
## [19] compiler_3.6.3          httr_1.4.2              backports_1.2.1        
## [22] assertthat_0.2.1        Matrix_1.3-2            lazyeval_0.2.2         
## [25] cli_2.2.0               tweenr_1.0.1            htmltools_0.5.1        
## [28] prettyunits_1.1.1       tools_3.6.3             gtable_0.3.0           
## [31] glue_1.4.2              rnaturalearthdata_0.1.0 maps_3.3.0             
## [34] Rcpp_1.0.5              cellranger_1.1.0        vctrs_0.3.6            
## [37] transformr_0.1.2.9000   iterators_1.0.13        crosstalk_1.1.1        
## [40] xfun_0.20               rvest_0.3.6             lpSolve_5.6.15         
## [43] lifecycle_0.2.0         scales_1.1.1            hms_1.0.0              
## [46] yaml_2.2.1              gridExtra_2.3           stringi_1.5.3          
## [49] highr_0.8               e1071_1.7-4             rlang_0.4.10           
## [52] pkgconfig_2.0.3         evaluate_0.14           lattice_0.20-41        
## [55] htmlwidgets_1.5.3       labeling_0.4.2          tidyselect_1.1.0       
## [58] plyr_1.8.6              magrittr_2.0.1          R6_2.5.0               
## [61] magick_2.5.2            generics_0.1.0          DBI_1.1.0              
## [64] pillar_1.4.7            haven_2.3.1             foreign_0.8-75         
## [67] withr_2.3.0             units_0.6-7             modelr_0.1.8           
## [70] crayon_1.3.4            KernSmooth_2.23-18      rmarkdown_2.6          
## [73] progress_1.2.2          readxl_1.3.1            data.table_1.13.6      
## [76] reprex_0.3.0            digest_0.6.27           classInt_0.4-3         
## [79] webshot_0.5.2           R.cache_0.14.0          R.utils_2.10.1         
## [82] munsell_0.5.0